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Abstract 

Factorial moments are convenient tools in particle physics to characterize the multi- 
plicity distributions when phase-space resolution (A) becomes small. They include all 
correlations within the system of particles and represent integral characteristics of any cor- 
relation between these particles. In this letter, we show a direct comparison between high 
energy physics and quantitative finance results. Both for physics and finance, we illustrate 
that correlations between particles lead to a broadening of the multiplicity distribution and 
to dynamical fluctuations when the resolution becomes small enough. From the generating 
function of factorial moments, we make a prediction on the gap probability for sequences 
of returns of positive or negative signs. The gap is defined as the number of consecutive 
positive returns after a negative return, thus this is a gap in negative return. Inversely for a 
gap in positive return. Then, the gap probability is shown to be exponentially suppressed 
within the gap size. We confirm this prediction with data. 



1 Introduction 



Multiplicity distributions and correlations between final- state particles in nuclear interactions 
are an important testing ground for analytic perturbative theory, as well as for Monte-Carlo 
(MC) models describing the hadronic final state [1J. Two-particle angular correlations have 
been extensively studied experimentally [2]. Specific statistical tools, namely the normalized 
factorial moments, have emerged in order to analyze in much details multiplicity distributions 
measured in restricted phase-space regions. The normalized factorial moments are defined as 

F q (A) = (n(n -l)...( n -q + l))/(n) 9 , q = 2, 3, . . . , 

for a specified phase-space region of size A. The number, n, of particles is measured inside A 
and angled brackets (...) denote the average over all events. The factorial moments, along with 
cumulants [3] and bunching parameters 01, are convenient tools to characterize the multiplicity 
distributions when A becomes small. Indeed, for uncorrelated particle production within A, 
Poisson statistics holds and F q — 1 for all q. Correlations between particles lead to a broaden- 
ing of the multiplicity distribution and to dynamical fluctuations. In this case, the normalized 
factorial moments increase with decreasing A. This effect is frequently called intermittency Q. 
As a matter of fact, it has been noticed in [5] that the use of factorial moments allows to extract 
the dynamical signal from the Poisson noise in the analysis of the multiplicity signal in high 
energy reactions. 

In addition, it has been shown that it is possible to define and compute a multi-fractal di- 
mension, T> q , for the theory of strong interactions (6l|7]| 

F q (A) = (n(n - 1) ... (n - q + 1)) / (n) q oc A {q - lKl -^ /d) (1) 

where d is the dimension of the phase space under consideration (d = 2 for the whole angular 
phase space, and d — 1 if one has integrated over, say the azimuthal angle). In the constant 
coupling case V q is well defined and reads 

F> q = 7o (2) 

q 

where 7q = ACaO!s/2it, as is the strong interaction coupling constant, Ca is the gluon color 
factor B6l|71|. The choice of the factorial moments as a specific tool in order to study the scaling 
behavior of the high energy multiplicity distributions is then useful to analyze the underlying 
dynamics of the processes. In principle, we can extend this last idea to other fields where 
factorial moments can be defined. 



2 Generating Function for Factorial Moments 

The multiplicity distribution is defined as P n = cr n /Y^L a ™' wnere a n is the cross section of 
an n-particle production process (the so-called topological cross section) and the sum is over 
all possible values of n so that 

oo 
n=0 



1 



The generating function can be defined as 

G(z) = £p B (l + z) B , 



(4) 



n=0 



which substitutes an analytic function of z in place of the set of numbers P n . Then, we obtain 
the factorial moment or order q as 



1 d"G(z) 



q (n)i dxfl 
and the corresponding definition for cumulants 



2=0, 



K„ 



1 dHnG(z) 



(n)« 



dz'i 



Z=0, 



(5) 



(6) 



The expression for G(z) can then be re-written as 



q=0 q - 



(7) 



lnG(z) = £fUn>% (#i = l). 

9=1 ? ' 



(8) 



The physical meaning of these moments has been discussed in the previous section. Another 
interpretation can be seen from the above definitions if they are presented in the form of integrals 
of correlation functions. Let the single symbol y represent all kinematic variables needed to 
specify the position of each particle in the phase space volume A [8]. A sequence of inclusive 
g-particle differential cross sections d g a/dy 1 . . . dy q defines the factorial moments as 



F q ~ 



dyi 



d q a 

a dy x ... dy q 



(9) 



Therefore, factorial moments include all correlations within the system of particles under con- 
sideration. They represent integral characteristics of any correlation between the particles 
whereas cumulants of rank q represent genuine g-particle correlations not reducible to the prod- 
uct of lower order correlations. 



3 Experimental Applications 

In Ref. J9l, it has been shown that factorial moments can be applied conveniently to quantitative 
finance. Namely, if we divide price series y(t) in consecutive time windows of lengths A, we 
can define a set of events. In each window, we have a certain number of positive returns n + , 
where y(t) — y(t — 1) > 0, and similarly of negative returns n_. If the sequence of returns is 
purely uncorrelated, following a Gaussian statistics at all scales, we expect F q — 1 for all q. 
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As explained in previous sections, correlations between returns may lead to a broadening 
of the multiplicity distributions (n + or n_ or even a combination of both) and to dynamical 
fluctuations. In this case, the factorial moments may increase with decreasing A, or increasing 
the number of bins that divide A. In Ref. flU, we have used the factorial moment of second 
order F 2 for like-sign returns 

2 ii «»+>/»-) 2 (10) 

where (n) is the average number of positive returns in the full time window (A), nu ns denotes 
the number of bins in this window and is the number of positive returns in k-th bin. Similarly, 
we have defined the unlike- sign returns for F 2 

p+- __ 1 J2k=T { n k U k } / n bins _j 

2 N ™ents J^ ts {n + )(n-)/{n bins ) 2 

with similar notations. Considering different price series, we have shown that for a small res- 
olution in time window 0, below 4 to 8 hours, a deviation with respect to pure Gaussian 
fluctuations is observed from the shape of second order factorial moments (flOl) and CCD) 0. 
Results are reminded in Fig. [TJ 

In order to compare the strength of the intermittency obtained in financial data, it is interest- 
ing to provide a comparison with standard high energy phenomenon in particle physics [|8l[T0ll. 
In what follows we discuss the results obtained with two of the most widely used e + e~ MC 
event generators, JETSET [[TT| and HERWIG |fi~2l. As it is done in the finance case, we define 
moments for like-charge and unlike-charge combinations of particles separately. Fig. |2] shows 
the behavior of F 2 ++ ' + (nbi ns ) at a centre-of-mass energy 91.5 GeV ifTOll . Without entering into 
details, we observe the rise of both observables with n bins = b as in Fig. [Q 

By comparing Fig. [Hand Fig. [2l we notice also that the intermittency is weak in the financial 
data w.r.t. high energy data, which is a reasonable observation. 



4 Gap Probability 

We can use the formalism introduced in section |2] in order to derive some further statements. 
From Eq. Q, we get 

G(-l)=po 

which corresponds to the probability to have zero particle in a phase space A or to have zero 
positive (resp. negative) returns in a given time window (finance case). This defines a gap, 
either in rapidity for particles or in duration for positive (resp. negative) returns. Let us use Eq. 
© to express G(—l) in another way using cumulants K q 

G(-l) =exp(-K 1 + K 2 /2\ + ...) 

When we can neglect correlations within a large time window flU, we have K 2 = and we 
conclude 

p ~ exp(— K\) = exp(— < n +i _ >) = exp(— pA) (12) 
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where A represents either the phase space or the time window. This last expression is very 
simple and instructive. It states that the gap probability is exponentially suppressed either in 
rapidity or in time. This is illustrated in Fig. [3]for the finance case. We observe the distribution 
of events (probability distribution) as a function of the size of the gap. The gap is defined as 
the number of consecutive positive returns after a negative return, thus this is a gap in negative 
return. Inversely for a gap in positive return. This gap is given in number of time units for the 
financial time series considered. In Fig. [31 we display results for the Euro future contract series 
over 10 years, sampled in two different time units. In both cases we observe effectively an 
exponential fall of the probability distribution as a function of the gap size. As illustrated also 
in Fig. [3l this exponential fall does not depend on the sampling. This confirms the prediction of 
Eq. (fT2l) . The case of particle physics is more complex, see Ref. lfT3l . 

5 Conclusion 

In this letter, we have discussed a direct comparison between high energy physics and quanti- 
tative finance results on factorial moments analysis. Both for physics and finance, we have il- 
lustrated that correlations between particles lead to a broadening of the multiplicity distribution 
and to dynamical fluctuations when the resolution becomes small enough. From the generating 
function of factorial moments, we have shown that p ~ exp(— pA). This expression states 
that the gap probability is exponentially suppressed within the gap size. The gap is defined as 
the number of consecutive positive returns after a negative return, thus this is a gap in negative 
return. Inversely for a gap in positive return. We have confirmed this prediction with data. 
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Figure 1: Euro future contract. Like-sign and unlike-sign factorial moments F 2 ++ and F 2 + ~ 
are displayed for 1, 2, 4, 10 or 20 bins. The global time window of 16.5 hours {nu ns = 1). 
This provides a time resolution ranging from 1.6 hours for 20 bins till 16.5 hours for 1 bin. We 
observe that for 1 and 2 bins segmentation, F 2 ++ is found equal to 1 and F 2 ++ is increasing above 
1 when the number of bins gets larger than 4. This confirms that non-Gaussian fluctuations in 
the sequence of returns returns start to play a role when the resolution is below 4 hours. 




b b 
o) b) 

Figure 2: e + e~ collisions at high energy. Factorial moments F 2 (y, <p) for like-charge (a) and 
unlike-charge (b) factorial moments F 2 ++ and -F 2 + ~, calculated using JETSET 7.4 for a center 
of mass energy of 91.5 GeV. See text for details. 
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Figure 3: Left: Euro future contract sampled in 5min time units. Right: Euro future contract 
sampled in 30min time units. The distribution of events (probability distribution) is displayed 
as a function of the size of the gap. The gap is defined as the number of consecutive positive 
returns after a negative return, thus this is a gap in negative return. Inversely for a gap in 
positive return. The Euro future contract series over 10 years is sampled in two different time 
units (5min-Left and 30min-Right). In both cases, we observe effectively an exponential fall 
of the probability distribution as a function of the gap size. This confirms the relation derived 
in this letter p ~ exp(— K{) = exp(— < n+ >) = exp(— pA). We have presented an 
exponential fit on top of each plot. On the Right plot (30 minutes sampling), we super-impose 
also the fit obtained on the 5 minutes sampling, in order to show that the exponential slope is 
identical in both cases. 
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